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Weighted  L  Quantile  Distance  Estimators 
For  Randomly  Censored  Data 

by 

R.  L.  Eubank  and  V.  N.  LaRlccia 

Southern  Methodist  University  and  University  of  Nebraska-Lincoln 

Short  Title:  Estimators  for  Censored  Data 

Summary:  The  asymptotic  properties  of  a  family  of  minimum  quantile 

function  distance  estimators  for  randomly  censored  data  sets  are  con¬ 
sidered.  These  procedures  produce  an  estimator  of  the  parameter 

2 

vector  that  minimizes  a  weighted  L .  distance  measure  between  the 
Kaplan-Meier  quantile  function  and  an  assumed  parametric  family  of 
quantile  functions.  Regularity  conditions  are  provided  which  insure 
that  these  estimators  are  consistent  and  asymptotically  normal.  An 
optimal  weight  function  is  derived  for  single  parameter  families, 
which,  for  location/scale  families,  results  in  censored  sample  analogs 
of  estimators  such  as  those  suggested  by  Parzen  (1979a,  1979b),  and 
Weiss  and  Wolfowitz  (1970).^. 

1.  Introduction.  In  this  paper  we  consider  the  problem  of 

parameter  estimation  from  randomly  censored  data  sets.  A  general 

method  of  estimation  is  presented  for  cases  when  the  data  is  assumed 

to  be  from  a  known  parametric  family.  The  technique  is  based  on  the 

2 

minimization  of  a  weighted  L  distance  measure  between  the  Kaplan- 

Meier  empirical  quantile  function  and  the  assumed  parametric  family 

of  quantile  functions  and  is  applicable  to  most  common  distributions. 

Let  X. denote  the  true  survival  times  of  n  individuals 
I  n 


2 


which  are  assumed  to  be  a  random  sample  from  the  distribution  function 
(d.f.)  F(x;j}°)  where  _0°  is  a  fixed,  possibly  unknown  element  of  a  known 
set  or  region  ©Ck  .  Further  let  Y^,...,Yn  denote  n  independent 
identically  distributed  censoring  random  variables  with  common  distri¬ 
bution  function  H  that  are  also  assumed  to  be  independent  of  the  X^'s. 
In  the  random  censoring  model  one  observes  not  the  X^'s  but,  instead, 
the  pairs  of  random  variables  where  Z^  =  min  (X^.Y^)  and 

5^  ■  Ijx  <  y  }  with  I  denoting  the  indicator  function.  The  d.f. 
of  the  Z^'s,  F*, is  then  given  by  the  relation 

(1.1)  1-F*  (x;9°)  -  [1-F(x;6°)]  [l-H(x)  ] . 

An  important  problem  associated  with  this  model  is  the  estimation  of 
the  parameter  vector,  9°,  from  the  observed  data. 

Let  F^(x)  denote  the  Kaplan-Meier  estimator  of  the  d.f.  F(x;9°) 
(Kaplan  and  Meier  (1958))  with  associated  empirical  quantile  function 
defined  by 

(1.2)  Q  (u)  -  inf {F  (x)  >  u}  . 

n  n  — 

x 

Estimation  problems  pertaining  to  specific  forms  for  the  vector  9^°  have 

been  addressed  by  Sander  (1975a,  1975b,  1975c)  and  Susarla  and  Van 

Ryzin  (1980)  and  Reid  (1981)  using  estimators  based,  explicity,  on 

both  F  and  Q  .  In  contrast,  we  develop  an  estimation  procedure 
n  n 

applicable  to  general  jJ°  under  the  assumption  that  the  functional 
form  of  F  is  known.  Thus  when,  for  instance,  9°  consists  of  only 
a  location  and  scale  parameter,  the  model  we  assume  is  the  censored 
sample  analog  of  the  classical  location  and  scale  parameter  model. 


Minimum  distance  estimation  procedures  based  on  the  d.f.,  the 
probability  density  function,  the  characteristic  function,  and  the 
quantile  function  have  been  proposed  and  both  their  large  and  small 
sample  properties,  for  non-randomly  censored  data  sets,  have  been 
extensively  investigated.  Specific  minimum  distance  procedures 
have  been  shown  to  possess  excellent  robustness  properties,  to  be 
consistent,  asymptotically  normal,  and  in  some  cases  fully  efficient. 
While  a  detailed  discussion  of  these  points  is  beyond  the  scope  of 
this  paper,  the  interested  reader  is  referred  to  Beran  (1977),  Parr 
and  Schucany  (1980),  and  Millar  (1981)  for  discussions  of  techniques 
which  utilize  the  d.f..  Estimation  procedures  formulated  in  the 
quantile  function  domain  may  be  found,  for  the  case  of  location 
and/or  scale  parameter  estimation,  in  Parzen  (1979a)  and  Eubank  (1981) 
and,  for  more  general  parametrizations,  in  LaRiccla  and  Wehrly  (1981) 
and  LaRlccia  (1981).  We  note  that  certain  results  presented  in  sub¬ 
sequent  sections  may  be  regarded  as  censored  sample  analogs  of  tech¬ 
niques  developed  in  Parzen  (1979a)  and  LaRiccia  (1981)  and  are  obtained 
through  the  use  of  work  by  Sander  (1975a,  1975b)  and  Reid  (1981)  on 
the  convergence  of  the  empirical  quantile  function  and  linear  functions 
of  order  statistics  for  randomly  censored  data. 

In  Section  2  we  define  our  estimator  and  present  our  principal 
results  regarding  its  asymptotic  properties.  The  proofs  are  provided 
in  Section  3.  Finally,  in  Section  4,  we  discuss  the  estimator’s 
efficiency  and  robustness  properties  and  provide  an  optimal  weight 
function  for  single  parameter  families. 

2 .  A  family  of  minimum  quantile  function  distance  estimators. 

For  each  6_  t  0  define  the  quantile  function  associated  with  F(x;9) 


(2.1)  Q(u;6)  -  inf{x:F(x;0)  _>  u}  . 

For  a  given  weight  function  W(u;0)  mapping  (0,1)  *8  +  1  and  associated 

real  valued  functions  W„  (0)  and  W„  (0),the  minimum  quantile  function 

S1  “  B2  " 

distance  estimator  of  the  true  parameter,  £  ,  is  any  vector  which 
minimizes 

02 

(2.2)  R(Q  ,0)  -  /  W(u; 0) [Q  (u)  -  Q(u;0)]2du 

n  —  n 


+  Wg  (DlQnttj)  ~  Q(B1;£)1  +  W02(i)tQn(62)  ‘  Q(02;-)J 

over  all  0^  e  0  where  0^  and  62  are  fixed  real  numbers  satisfying 
0  £  8^  <  0£  —  1  ^ut  are  otherwise  arbitrary.  Since  individual  members 
of  this  family  of  estimators  are  distinguished  by  their  specific  weight 

A 

function  we  adopt  the  notation  £(Qn,W)  for  this  estimator.  For 

notatlonal  convenience  it  will  be  useful  to  have  an  expression  which 

incorporates  all  three  of  the  functions  Wr  W  and  W  .  Therefore, 

B1  B2 

for  any  function  Z(u;0)  we  define 

(2.3)  /62W*(u;0)Z(u;0.)du  -  /SlW(u;£)Z(u;0)du 


+  J  wg  (0)z(si;0) 

i-l  1 

and  (2.2)  may  now  be  written  as 


(2.4)  R(Qn;0)  -  /  2  W*(u;0){Qn(u)  -  Q(u;0)]Zdu 

81 


Weight  function  selection  strategies  will  be  discussed  in  Section  4. 
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In  general,  the  computation  of  jKQ^W)  can  be  accomplished  using 
standard  iterative  techniques.  Estimator  computation  is  particularly 
simple  for  location  and  scale  parameter  families  when  the  weight 
function  is  independent  of  £  as,  in  this  instance,  the  estimator 
is  readily  seen  to  have  a  closed  form. 

One  advantage  of  the  minimum  distance  estimation  technique  is 

that  R(Q  ,8(Q  ,W))  provides  a  measure  of  the  goodness-of-fit  of  the 
n  —  n 

assumed  parametric  family.  In  addition,  the  quantile  function 
approach  has  the  consequence  that  the  estimation  procedure  is 
directly  related  to  Q-Q  plotting  techniques.  Therefore,  the  esti¬ 
mators  can  be  easily  incorporated  into  a  statistical  package  which 
not  only  estimates  the  parameters  but  also  provides  checks  for  the 
appropriateness  of  the  assumed  parametric  family  and  graphs  which 
can  be  employed  to,  perhaps,  suggest  a  more  suitable  family  of 
distributions. 

A 

The  asymptotic  behaviour  and  distribution  theory  for  jKQn»W) 
is  the  subject  of  Theorem  1.  To  prove  these  results  we  require 
certain  restrictions  on  both  W  and  Q.  Therefore,  let  G  denote  the 
class  of  left  continuous  functions  on  (0,1)  that  are  of  bounded 
variation  on  (y,1-y)  for  all  0  <  y  <  and,  for  h,  g  e  G  define 


d(g,h)  -  sup|h(x)-g(x) |  . 

Also  for  any  function  Z(u;6)  let  Z*(u;8)  *  —  , 

2  ~  '  *  ‘ 

Z13(u;e)  .  ■  ,  Z'(uj9)  -  and  Z.  *  2(-;S)  . 


ae13ej 


The  following  assumptions  are  then  required  for  Theorem  1: 


(1)  d(Qn,Qg0)  -»p0,  where  -*-p  denotes  convergence  in 
probability. 

(2)  Sn  d(Qn,Qg0)  »  Op(l) ,  where  0p  denotes  probability  order. 

(3)  Let  V  denote  the  vector  having  ith  component 

— n 

g 

2W*(u;i°)Qi(u;e°)(Q  (u)-Q(u;0°))du) 

S1 

where  the  integral  is  defined  as  in  (2.2).  Then 

converges  in  law  to  the  Ng(0,A)  distribution,  denoted 

V  -*-,N  (0,A),  where  the  ijth  element  of  A  is 
— n  L  s  — 

-Si  -  0  Q*(u;jJ°)0^  (v;6°) 


fQ(u;0°) fQ(v;0°) 


—  K(u,v)dudv 


(2.5)  a  -  /  2f  2  W*(u;e0)W*(v;6°) 

J  81  h 

with 

min(u.v)  ,  o-l 

(2.6)  K(u,v)  -  (l-u)(l-v)/  t(l-s)Z(l-HQ(s;0°))]  Xds  , 


HQ(u;0°)  -  H(Q(u;00)),  and  fQ(u;0°)  -  F’(Q(u;0°);  0°) 
is  the  density-quantile  function  corresponding  to  F  (c.f. 
Parzen  (1979a)) . 

(4)  For  fixed  0_,  the  functions  W(u;£),  W*(u;£),  VT^(u;0),  Q*(u;6) 

and  Qi;^<u;0)  are  continuous  on  and»  ^or  **xec*  u,  in 

a  neighborhood  of  0°.  W*  (0),  W^(0),  k  ■  1,2,  are  continuous 

Sk  \ 

in  0.. 

(5)  The  functions  W*(u;0),  W(u;0}Q*(u;0)  ,  W^(u;0),  Wi(u;6^), 
Q^(u;0)  and  W(u;0_)Q^  (u;0)  are  bounded  by  integrable  functions 
uniformly  for  all  _0  in  the  neighborhood  of  j3° . 


Theorem  1  Under  the  regularity  conditions  (1)  -  (6) 

(i)  As  n-x»  there  exists,  with  probability  tending  to  one,  a 
unique  function,  jHQn,W),  which  locally  minimizes  (2.1), 

A 

(ii)  9 (Q  ,W)  is  a  consistent  estimator  of  9°, 

—  n  — 

(iii)  vG(0(Q  ,W)  -  9°)  -*■  N  (0,C)  where 
—  n  —  s 

(2.8)  C  =  [B(Q0O,0°)]_1A[B(Q9O,0O)]“1 

with  A  and  B  as  defined  in  (2.5)  and  (2.7). 

Conditions  (1)  -  (6)  can  be  replaced  by  many  different  sets  of 
restrictions.  In  particular  conditions  which  imply  (1)  and  (2)  can 
be  derived  from  the  results  of  Sander  (1975a)  whereas  restrictions 
implying  (3)  can  be  found  in  Sander  (1975b)  and  Reid  (1981).  As  an 
illustration  of  alternative  conditions  we  present  the  following 
corollary. 

Corollary  .  Let  assumptions  (4)  -  (6)  be  satisfied  and  further  assume 
that 

(i)  0  <  8^  <  <  1  and  H  is  a  continuous  function  satisfying 

HQ(62;l°)  <  1, 

(ii)  F  is  a  strictly  increasing  function  with  F(O;0_°)  *  0  , 

(iii)  W(u;0^°)Qi(u;0°)  is  differentiable. 


Then,  the  conclusions  of  Theorem  1  hold. 
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Before  proving  Theorem  1  in  the  next  section  we  note  that  an 
alternative  approach  to  our  method  of  proof  could  be  taken  wherein 

A 

8 (Q  ,W)  is  treated  as  a  functional  on  the  metric  space  (G,d).  Thus, 

—  n 

assumptions  similar  to  those  in  Beran  (1977)  on  the  differentiability 
of  the  functional  jKQQ*W)  could  be  used  to  replace  conditions  (4)-(6) 
although  the  present  conditions  seem  easier  to  verify  in  practice. 

However,  if  W(u;0°)  >  0  for  all  u  e  It  Is  usually  easier  to 

g 

check  the  condition  Z  a.Q*(u;6°)  j  0  for  all  a  *  (a.,..., a  )C  i*  0 
i-1  1  -Is 

which,  in  this  case,  is  equivalent  to  (6).  We  note  that  this  latter 

condition  is  an  assumption  on  the  identlf lability  of  the  parameters. 

3R(Qn,6) 

3.  Proofs.  We  now  prove  Theorem  1.  Let  D(Q  ,0)  *  - rr - 

- .  n  —  a v 

a  R«L.!>  th 

and  J(Q  ,0)  ■  - s -  .  Then,  D  is  the  s  x  1  vector  with  i 

n“  302 

element 

ft 

-  / ^2{(QnCu)-Q(u;0))2W*i(U;0)  -2W*(u;0)Qi(u;0)(Qn(u)-Q(u;0))}du 

S1 

*i 

where  the  W  notation  and  the  integrals  are  defined  in  a  manner  analogous 

to  (2.2).  Similarly,  J  is  the  s  x  s  Jacobian  matrix  with  ijth  element 

Rij(Q  ;0)  -  /  2{(Q  (u)-Q(u;0))2W*ij(u;0)-2(Qn(u)-Q(u;0))[W*i(u;0)Qd(u;£) 
n  —  in  —  —  n  — 

P1 

+  W*d(u;£)Qi(u;£)4W*(u;0)QiJ(u;9.)]+2W*(u;0)Qi(u;£)Qd(u;£)}du  . 

Now,  employing  the  compactness  of  and  conditions  (4)  -  (5)  tedious 

but  straightforward  calculations  give  the  following  three  results: 

(A)  In  an  open  neighborhood  of  (QgO,£°) ,  R^(Qn>£)  and  R^’d(Qn*j0) 
are  continuous  functions  with  respect  to  the  metric  p,  defined 
for  g,  heG  and  £,ae@  by 

p[(h,0),(g,a)]  ■  max{d(g,h) ,  I  e1"ot1l » *  *  * » 1 0s"a8 1  > 


9 
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(B)  /n/  2(0  (u)-Q(u;J»0)) V^CujJ^Jdu  ->-  0  and 

61  ! 

(C)  Rlj(Q  ,0)  -V  RiJ(Q  0,0°)  if  0  +  6°. 

n  —  o_  —  —  r  — 

Since  D(QgO,0°)  =  0  and  J(Qqo,0°)  is  positive  definite  it  follows, 
using  (A),  that  the  conditions  of  the  implicit  function  theorem  are 
satisfied.  Consequently,  there  exists  an  open  neighborhood  of  (Q  o,0°) 
and  a  unique  continuous  mapping,  6_(Qn»W) ,  of  an  open  ball,  8,  containing 

S  A  A  o 

Q  0,  into  R  such  that  d(Q  ,  0(Q  ,W))  *  0  and  9(Q_o,W)  *  0°.  Further,  as 
o_  n  —  n  —  — 

a  result  of  (1),  with  probability  tending  to  one,  Q  is  in  B  so  that 

n 

there  exists,  also  with  probability  tending  to  one,  a  unique  point 

A  A  A 

0(Q  »W)  such  that  D(Q  ,0(Q  ,W))«  0.  To  see  that  0(Q  ,W)  minimizes  R(Q  ,0) 
—  n  n -  n  — n  n— 

A 

note  that,  using  (A),  J(Qn»6.(Qn,w) )  ^p^^e0’  0.  )  which  is  positive 
definite  by  assumption  (6).  Thus,  result  (i)  has  been  shown.  Result 
(ii),  the  consistency  of  0(Q^,W),  is  an  immediate  consequence  of  the 

A  A 

continuity  of  0(*,W)  and  the  fact  that  0(Q0o,W)  *  _0°. 

By  the  mean  value  theorem 

D(Q  ,0°)  -  D(Q  ,0(Q  ,W))  +  J(Q  ,?)(0°-  0(Q  .W)) 
n  —  n  —  n  n  —  —  n 

where  0  -  0°  +  A(0°-0(Q  ,W))  for  some  sxs  diagonal  matrix.  A, 

—  —  —  —  n 

A 

with  all  elements  between  0  and  1.  Since  D(Qn».9(Qn»w))  ■  0, 

v£[i<Qn,W)-0°]  -  -J(Qn,0)'1^b(Qti,0°)  . 

Result  (iii)  now  follows  using  (B) ,  (C),  assumption  (3)  and  a 
Slutsky  type  argument. 

To  verify  the  corollary  we  first  note  that  it  is  shown  in 
Sander  (1975a,  Corollary  1)  that  */n[Qn(u)-Q(u;0o)  ]  converges  weakly 
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to  (-l/fQ(u;0))X(u)  on  where  X(*)  is  a  zero  mean  Gaussian 

process  with  covariance  kernel  (2.5).  Thus,  it  can  be  readily 

shown  that  assumptions  (1)  and  (2)  are  satisfied.  This  fact  also 

implies  that,  for  i=l,...,s, 

2 

I  w  (0o)Qi(8k;0°)(Q  (e.)  -  Q(ek;0°))+  N(o,a*) 

k*l  °k  LI 

where 


Since  W(u;0^°)Qi(u;0^°)  is  differentiable  on  the  results  of 

Example  3  of  Reid  (1981)  have  the  consequence  that 

fj  _ 

/  w(u;0°)Q  (u;0°)  (Q  (u)-Q(u;0°))du  N(O,0_) 

S1  n  "  L  2 

where 

8  8 

a2  “  /  2/  ^(v;£°)Qi(v;£°)W(u;£°)Qi(u;e^c>)  [fQ(v;e^°)fQ(u;^°)]_1X(v,u)dvdu  . 
B1  B1 

Then,  an  argument  similar  to  that  used  in  Shorack  (1972)  gives  (3)  and 
the  conditions  of  Theorem  (2.1)  are  satisfied. 

4.  Robustness  and  efficiency.  In  this  section  we  discuss 

the  robustness  and  efficiency  properties  of  the  estimator  0(Q  ,W) . 

—  n 

Using  an  approach  similar  to  that  of  Parzen  (1979a)  an  optimal  weight 
function  is  provided  for  estimation  in  single  parameter  families. 

In  the  case  of  location  or  scale  families  with  no  censoring  (H»0)  the 
resulting  estimators  agree  with  estimators  proposed  by  Parzen  (1979a,  1979b) 
and,  consequently,  are  similar  to  those  given  by  Weiss  and  Wolfowltz  (1970). 
The  results  in  this  section,  therefore,  provide  an  extension  to  randomly 
censored  samples  of  the  work  of  these  authors. 

It  follows  from  the  proof  of  Theorem  1  that,  asymptotically, 

A 

0(Qn»W)  is  a  weighted  sum  of  s  linear  functions  of  the  order  statistics. 
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Thus  the  Influence  curve  for  these  estimators  can  be  determined  using 
the  techniques  developed  in  Reid  (1981).  For  a  given  parametric  family 
of  distributions,  weight  functions  with  specific  properties  or  which 
provide  specific  types  of  protection  can  then  be  derived.  In  parti¬ 
cular,  by  proper  selection  of  W(u;j3)  an  estimator  which  is  locally 
robust  or  one  that  has  specified  efficiency  properties  can  be  obtained. 

The  remainder  of  this  section  is  devoted  to  the  problem  of  weight 
function  selection.  To  obtain  an  optimal  weight  function  we  first 
note,  as  in  Parzen  (1979a),  that  the  covariance  kernel,  K,  in  (2.6) 
generates  a  reproducing  kernel  Hilbert  space,  ff(K) ,  which  provides  us 
with  a  natural  measure  of  information  on  The  properties  of 

H( K)  that  impact  on  our  present  objective  will  now  be  briefly  dis¬ 
cussed.  The  reader  is  referred  to  Aronszajn  (1950)  and  Parzen  (1961a, 
1961b)  for  a  more  detailed  presentation  of,  respectively,  the  theory 
of  reproducing  kernels  and  their  role  in  inference  for  stochastic 
processes. 

Assuming  that  H  admits  a  continuous  density,  h,  it  follows  from 
Sacks  and  Ylvisaker  (1966)  that  H( K)  consists  of  continuous  functions 
on  {8^,62!  having  finite  H( K)  norm  where,  for  geff(K),  the  norm  is 
given  by 

g 

llgll*  *  /  2[g'(u)]2[l-HQ(u;_9)]du 
K  B, 


(4.1) 


2 

f2  P g(u)1  hQ(u;£)  , 

{  (  )  fQ<«;ir du 

, r  1  1-HQ(B- ie) 

+  g(8  )2 - — 

1  I  K(B1,B1)  1  -  Bx 

g(s2)2 


+ 


U-hq(b2;!)] 
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The  Inner  product  in  H( K)  will  be  denoted  by  <•»•>£.  When  considered 
as  a  function  of  s  (for  fixed  t)K(s,t)  possess  the  so  called  reproducing 
property  that  for  geff(K) 

(4.2)  <g,  K(‘,t)>K  -  g(t). 

We  now,  and  in  subsequent  discussions,  Impose  the  regularity 
condition  that  fQ(u;.0)Q*(u;0.)  e  H( K)  and  define  the  information 
measure 

(4.3)  1(0^  -  ||fQ(-;0)Q1(*;0)||*  . 

To  justify  this  definition  we  first  note  that  the  densities  corre¬ 
sponding  to  the  censored  and  uncensored  observations  are,  respectively, 
fC(x;0)  -  [1-F(x;0)]h(x) 

and 

fu(x;0)  *  [1-H(x)]f(x;£)  . 


Then  by  making  the  change  of  variable  X  *  Q(u;6)  it  can  be  shown  that, 
for  instance,  when  8^  ■  0  and  &2  m  1*  becomes 


(4.4) 


KV  -  / 

o 


fU(x;J3)  dx 


+ 


fC(x;0)dx 


This  is  now  recognized  as  the  Fisher  information  corresponding  to  the 
parameter  0^  upon  examination  of  the  form  of  the  likelihood  for  the 
Z^'s  given  in,  for  example,  Kalbf leisch  and  Prentice  (1980).  Similarly, 
1(0^)  can  be  seen  to  provide  a  measure  of  information  in  the  Fisher 
sense  when  0  <  8^  <  8j  <  1. 
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Now  suppose  that  there  is  only  one  unknown  parameter,  0^  ■  6, 
and  let 


G(u;6)  *  (1-u)/  - = — -  du»  . 

0  (l-oi)  [1-HQ(oj;6)  ] 

Under  the  assumption  that  fQ(u;0)Q^(u;0)  is  twice  continuously 
differentiable  define 


(4.5)  4> (u; 0)  =  -  [fQ(u;0)Q1(u;0)  ]"[l-HQ(u; 0)  ] 


_  poQo(u;e)ql(u;9)  J  (1.u) 


hQ(u;0) 

fQ(u;6) 


for  8^  <  u  <  with 

{fQ(B1;0)Q1(81;0)G,(8.;0)  -  G(B. ;0) (fQ-Q1) ' (8, ;0) > [1-HQ(B, ;0) ] 

(4.6)  6  (0)  -  - 1 ± - 1 - 1 - 1 - ^ - 

B1  G(B1;0) 


and 


(4.7) 


*  6,(9) 


{fQ(B2;9)Q1(B2;0)  +  (l-B2)(fQ*Q1),(B2;0)}[l-HQ(B2;0)] 


where  (fQ-Q1) » (Bk;0)  -  [fQ(u;0)Q1(u;0) ] ' 

from  Sacks  and  Tlvisaker  (1966)  it  follows  that 
the  representation 


Then,  using  results 
fQ(u;0)Q1(u;0)  admits 


(4.8)  fQ(u;0)Q1(u;0) 


for  u  c  [B1,B21.  By  using  (4.8)  in  conjuction  with  the  reproducing 
property  for  K  it  is  possible  to  verify  the  two  important  identities 
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||fQ(-;0)Q1(-;e) 

(4.9) 


Z 

■  /  fQ(u;6)Q1(u;6)<j)(u;9)du 

61 

2 

+  l  (9)fQ(Bk;9)QX(Bk;9) 
k-1  k 


and 

S  8 

|  |fQ(*;9)Q1(*;9)  |  |^  =  /  2/  2<t>(u;9)<j>(v;9)K(u,v)dudv 

»isi 

(4.10)  2  B, 

+  2  1*  (9)/  X  *(u;0)K(B.  fu)du 

B,  k 


2  2 

+  M  ♦«  (e)*»  (0)K(Bv»eo) 

k-1  1=1  \  £ 


Now  define  the  weight  function 

(4.11)  »(u;8)  -  iiaiy.SKsi?)  , 

Q  (u;9) 


with 

(4.12) 


W  (9) 
6k 


^g^(9)fQ(8k;9) 

Q1(0k;e) 


Bx  <  u  <  B2 


k  -  1,2. 


Using  this  weight  function  in  (2.4)  and  (2.6)  and  applying  the 
identities  (4.9)  -  (4.10)  it  can  be  concluded,  from  Theorem  1,  that 


i^n(0(Q  ,W)  -  9°)  -*■  N(0, 1/1(9°)) 
n  L 

and  hence,  the  weight  function  (4.11)  -  (4.12)  is  optimal  for  the 
estimation  of  9°. 

The  presence  of  l-HQ(u;9)  in  the  weight  function  has  the  con¬ 
sequence  that  the  weight  function  depends  on  the  censoring  distribution 
which  may  not  be  known.  Although  the  details  are  beyond  the  scope 
of  this  paper,  such  difficulties  can  be  avoided  by  using  strongly 
consistent  estimators  of  H  and  h  in  W  and  W0  ,  k  -  1,2. 
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To  provide  further  Insight  into  the  properties  of  this  weight 
function  consider  the  special  case  when  9  is  a  scale  parameter  and 
there  is  no  censoring,  i.e.,  H  ■  0.  In  this  case,  we  have 

F(x;9)  -  F(x/0) 
f(x;9)  -  -|f(x/9) 

Q(u;6)  -  0Q(u) 
fQ(u;0)  -  |  fQ(u) 

and 

Q1(u; 0)  =  Q(u) 

where  F(x) ,  Q(u)  and  fQ(u)  are  all  specified  functions.  Equations 
(4.5)  -  (4.7)  now  become 


(4.13)  *(u)  -  -(fQ(u)Q(u) }" 

(4.14)  ♦  -  ~  fQ(81)Q(61)-(fQ-Q)’(B1) 

and 

(4.15)  ♦  -  fQ(B2)Q(B2)  +  (fQ-Q)’(B2) 


so  that  W  is  Independent  of  0.  Consequently,  0(Q  ,W)  may  be  expressed 

n 


in  the  closed  form 
(4.16)  0(Qn,W)  - 


S2  2 

L  *(u)fQ(u)Q  (u)du  +  r  fQ(6.)Q  (B,  ) 
P1  k-1  Bk  k  n  k 


/  2  <Ku)fQ(u)Q(u)du  +  l  4>r  fQ(B.)Q(B.) 
Bx  k-1  Bk  k 


This  is  precisely  the  fully  efficient  estimator  of  0  considered  by 
Parzen  (1979a, 1979b) .  Identical  results  hold  for  location  parameter 
estimation  as  well.  Parzen  has  noted  that  these  estimators  are 
similar  to  those  presented  in  Weiss  (1964)  and  Weiss  and  Wolfowitz  (1970) . 
Consequently,  the  minimum  quantile  function  distance  procedures  developed 
in  this  paper,  when  incorporated  with  the  weight  functions  derived  in  this 
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